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Abstract 

We perform a quantum mechanical analysis of the pendular cavity, using the positive-P rep- 
resentation, showing that the quantum state of the moving mirror, a macroscopic object, has 
noticeable effects on the dynamics. This system has previously been proposed as a candidate for 
the quantum-limited measurement of small displacements of the mirror due to radiation pressure, 
for the production of states with entanglement between the mirror and the field, and even for super- 
position states of the mirror. However, when we treat the oscillating mirror quantum mechanically, 
we find that it always oscillates, has no stationary steady-state, and exhibits uncertainties in po- 
sition and momentum which are typically larger than the mean values. This means that previous 
linearised fluctuation analyses which have been used to predict these highly quantum states are of 
limited use. We find that the achievable accuracy in measurement is far worse than the standard 
quantum limit due to thermal noise, which, for typical experimental parameters, is overwhelming 
even at 2 mK. 

PACS numbers: 42.50.Lc,42.50.-p,42.50.Pq,03.65.Yz 
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I. INTRODUCTION 



The pendular cavity, or Fabry- Perot cavity in which one of the mirrors is free to oscillate, 
has previously been investigated by a number of researchers, both experimentally 
and theoretically . Closely related schemes have been theo- 

retically proposed to entangle mirrors jl4] and create superposition states of a mirror isl ]. 
Common to almost all the theoretical treatments is a linearisation of quantum Langevin 
equations around their classical steady-state solutions Qj. This then allows for the rela- 
tively simple calculation of spectral quantities which may be measured outside the cavity. 
Results obtained in this way have been used to analyse the sensitivity of gravity wave in- 
terferometers [l7|, predict the suppression of quantum noise 0,0], propose quantum nonde- 
molition measurements of photon number |6|], analyse the quantum limits to measurements 
with an atomic force microscope analyse the quantum noise in position measurements 
of the oscillating mirror Q], calculate the phase noise in the cavity field predict the en- 
tangleme,. of .naooscop.c osciUaprs via radiation pressure Q. and propose the quantum 
locking of interferometer mirrors [13|| . Using the state- vector approach so common in quan- 
tum computing theory, it has been proposed that quantum superpositions, entanglement 
and near-number states of the cavity field, along with superposition states of the mirror, 
can be produced with this system ^8]. Using a similar state- vector approach, it has been 
proposed that quantum superpositions of a mirror may be created by the interaction with 
a single photon 

It is well-known that the linearised fluctuation analysis used in the majority of the the- 
oretical papers cited above is limited in its applicability. It has been shown, for example, 
that the mean-field equations derived in this way can give misleading results for travelling- 



wave second harmonic generation 



light and condensed atoms 
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and for the intracavity interaction between 



, .221. 12311 as well as in Raman photoassociation of atomic 

nrrTTi 

Bose- Einstein condensates |2^ |2^ |2a |22l| . The spectra calculated via this method are also 
known to not be accurate near any critical points of the system, as has been shown with 



the optical parametric oscillator [28|, |29| . There are two conditions which must be fulfilled 



for a linearised analysis to be trustworthy. The first has to do with the sign of the real 
part of the eigenvalues of the drift matrix of the equations written for the fluctuations; if 
these have the wrong sign the fluctuations can grow exponentially and the analysis loses its 
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validity. The second has to do with the size of the fluctuations themselves, in relation to the 
classical steady-state values. As we will show, this second condition is not fulfilled for this 
system, as the quantum state of the mirror, or more accurately, of the mirror phonons, is 
reasonably expected to be thermal. A characteristic of a thermal state is that the variance is 
larger than the mean value, which therefore makes any expression of the mirror phonons as 
haviner some well-defined classical mean value plus small fluctuations rather dubious. Even 
though cooHn. of th. ...o. v,a feedback .echanis,ns ha. bee„ achieved Q. a„d analysed 
theoretically j3l|, due to the nature of the coupling between the electromagnetic field and 
the mirror phonons, it seems that all that can be achieved is a thermal state at a lower 
temperature, so that the problem remains. 

To treat the macroscopic mirror quantum mechanically we will begin with the Hamilto- 
nian approach of Law js^, in the approximation that only a single optical mode is important, 
and extended to include cavity pumping and damping. To treat the fluctuations of the mir- 
ror which result from its coupling to a thermal reservoir, we will use the Brownian motion 

. n 

master equation developed by Diosi which is suitable for the temperatures we will 
consider here. Following a common procedure in quantum optics [34j, we will develop a 
Fokker-Planck equation in the positive-P representation [35j. This Fokker-Planck equation 
allows us to write stochastic differential equations which are an exact mapping from the 
system master equation, and which can be used to calculate any desired normally-ordered 
operator moments. 



II. SYSTEM AND HAMILTONIAN 



We consider a system of a pumped Fabry-Perot cavity in which one of the mirrors is free 
to undergo oscillatory motion due to both the light pressure and thermal fluctuations. We 
use the standard annihilation, a, and creation, operators for the electromagnetic field, and 
the operators x and p for the displacement from the equilibrium position and the momentum 
of the mirror, which will be treated as a harmonic oscillator. Neglecting the coupling of the 
mirror to its bath for the moment, we can write the Hamiltonian as 

H = Hfree + -f^int + Hp^mp + Hi^athi (l) 
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where 

1 1 

Hint = —hga^ax, 

Hbath = To) + V^a. (2) 

In the above, ujq is the field frequency, e represents the classical real pump, m is the mass 
of the mirror, uOm is the mirror oscillation frequency and g = ujq/ L is the coupling between 
the mirror and the cavity field, with L being the length of the cavity. The Fs represent 
optical bath operators. The damping of the mirror, which we will treat as Markovian, will 
be included at the next step. 

We now wish to write a master equation for the density matrix of our combined system in 
a frame rotating at ujq. To do this, we will make two different, but consistent approximations 
for the damping of the cavity and the mirror. The cavity reservoir will be considered to 
be at zero temperature, which is consistent with the very high temperatures necessary to 
produce thermal photons at the frequencies involved. The mirror reservoir will be treated as 
being at a finite temperature, which is necessary because of the number of thermal phonons 
which will be present in the system. As the temperatures required to create these respective 
excitations differ by many orders of magnitude, these approximations are not contradictory. 
This process gives us 



.,dp 



H,p +Cp 

- — p'^ + -mcu^x^ — hga)ax — ihe{d — a)), , 



+ih'j {2apa) — a)ap — pa^aj + -DmP, (3) 

where 7 represents the loss rate through the fixed mirror and D^p represents the mirror 
damping, using the Brownian motion master equation developed by Diosi js^. 



DmP = 7m [X, {p, p}] - [X, [X, p]] - ^^^'^^'^B [p, [p, p]] . (4) 

dB 



In the above, 7^ is the mirror damping rate, which depends on temperature through the 



mechanical quality factor, A^^ = h/y/AmkBT, the thermal de Broglie wavelength of the 
mirror, with /c^ Boltzmann's constant and T the temperature, k is a numerical factor 



which must be greater than 1 for this master equation to be of the Lindblad form, but is 
not important here as it leads to terms of the order hum/kBT, which we will show to be 
insignificant at the temperatures we consider. 



III. EXPANSION OF THE MIRROR IN COHERENT STATES 

Rather than writing the Heisenber g eq uations of motion, which are difficult to solve, we 
will make use of the original definition Isjl of the annihilation and creation operators in terms 

n 

of X and p and develop stochastic differential equations in the positive-P representation j35|| . 
This allows us to use c-number equations which describe all the quantum properties of the 
mirror dynamics contained in the original master equation. We will describe the operators 
X and p, in terms of the operators b and b\ where 



X = A(h + b^ 
p = B(h-py (5) 



with 

A 



2'muJri 



and [b, b^ = 1. Writing the equations using these variables is advantageous because it allows 
us to automatically define a P-representation of the density matrix in terms of an expansion 
in the minimum uncertainty (coherent) states |/?), defined as b\P) = It also means 

that the mirror quadrature variances have a coherent state or vacuum value of 1, the same 
as for the electromagnetic field. In fact, A and \B\ represent the standard quantum limits 
(SQL) for measurement of the mirror position and momentum, respectively. We note that 
these phonon annihilation and creation operators have previously been used to describe the 
mirror, but not in the context of developing phase-space representation stochastic differential 
equations 

In terms of these new variables, the master equation for the mirror damping is now 
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written as 



dB d2 



h 



b + b^ 



,P 



(7) 



IV. STOCHASTIC EQUATIONS 



Using the well-known operator correspondences for the P-representation [sUj, we may map 
the master equation onto a partial differential equation for the P-function of the system, 
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The diffusion matrix of the above equation is 
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(9) 



We note here that this drift matrix has diverging terms as T — 0, but this is not a problem as 
the Diosi master equation is valid in the limit where kBT ^ tiUm- As a physical example, in 
Ref. 0, wefindcj„ = 1.6x10^ s"^ so that hu^ = 1.72x10^29 whereas ksT = S.SxlO^^s j 
at 4.2 K, the temperature which we will mainly use in our investigations. 



If we wish to treat Eq. |S1 as a genuine Fokker-Planck equation which we may map onto 
stochastic differential equations, the matrix D must be positive-definite. Numerical inves- 
tigations using typical parameters show that this is not the case, therefore for quantum 
calculations we will have to use the positive-P representation 35j. The positive-P represen- 
tation equations in a doubled phase space can be found by the simple change of variables 
a* — i> a^, (3* — >■ (3^ , so that, (noting that a+ = a* only in the mean and similarly for 
13^)1 we now have four independent stochastic variables. Ignoring the terms proportional to 
hum/kBT due to their small relative size, one possible factorisation of the diffusion matrix, 
D = NN^, of Eq.iniis 
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which allows us to write a set of four stochastic differential equations (Note that the Ito 
form and the Stratonovich form of these equations are identical.). 
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In the above, the real Gaussian noise terms have the correlations 



(11) 



0, 



Vi{t)Vj{t') = 6i,6{t-t'). 



(12) 



The set of coupled equations may be integrated numerically, with averages taken over a 
large number of stochastic trajectories, which allows for the probabilistic calculation of any 



desired normally-ordered operator moments. As an example, with N trajectories, we have 

N 



a' -a-) = lim i-V" (13) 



where j labels the results from the jth trajectory. 

V. CLASSICAL ANALYSIS 
A. Steady-state solutions 

Before we return to the full stochastic equations, which we will solve numerically, we will 
investigate some of the classical properties of the system, which allow for analytical insights. 
From the drift part of Eq. |H1 we can immediately write the mean-field equations using the 
notation z for the classical mean-field value of z, 

_ = e--fa + tgAa{P + P*), 

dd _ _ 

^ = -tUJmf3 - 7n^(/5 - + igA\a\\ (14) 

from which we may find the classical steady state solutions. 

Solving Eqs. El for the steady-states, we find that (3ss is real, which means that the 
steady-state momentum is zero. (Note that this will not be the prediction of stochastic 
integration of the full equations.) However, using this fact we may write the solutions as 

R — R* — 1 2 

Pss Pss l^ssl ) 



e 



Although the solutions above are not closed (the solution for ass is a function of /J^s, etc.), 
we can make an iterative expansion, beginning with the result for fixed mirrors, 

oPss = (16) 
7 

and substitute this into the solution for f3ss- This can then be substituted into the solution 
for ass, the process being repeated until we attain the required degree of convergence. We 
note here that there are parameter regimes for which this expansion does not converge and 
that these are regions where we do not find classical steady-state solutions, but rather a 



limit-cycle, self-pulsing behaviour 



0. 
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Although we will demonstrate below that the classical steady-state solutions, especially 
for X and p, are not accurate in any parameter regime, they do allow for some insight into 
which property of the electromagnetic field outside the cavity is most likely to allow for 
an inference of the mirror position. The usual candidates are the intensity, which may be 
measured by photodetection, and the quadratures, which may be measured by homodyne 
detection. Defining the intracavity quadratures as Xa = a + a) and Ya = —i{a — a^), with 
their classical equivalents written in terms of Uss and q;*^, we find 

2e7 



" 2 I n'2rr2 ' 



7^ + 
2gex 



7^ + g^'x'^ 
7^ + g'^x'^ 

where we have set e as real and used x = 2A(3, with (3 real. As 7^ is typically much larger 
than g'^x^ , we may make a series expansion of these expressions. We find that 

7 V 1 

2^2 



2gex ( g X 



Ya^ ^\\ 



9 I 9 

7 V T 

It is immediately obvious that Ya depends on x to first order, while the other two exhibit only 
a second order dependence. This shows that homodyne measurements of the Ya quadrature 
will be more sensitive to variations in the position of the mirror than will be the other two 



measurements, as previously noted by Vitali et al. [3l|, although the Xa quadrature will 
show a weaker dependence rather than being totally independent of x as in the linearised 
analysis of Ref . jSl | . 

B. Bistability 

It has been predicted that, with a non-zero detuning between the field and the cavity res- 
onance, this system can exhibit bistability in the optical intensity Q]. To find the condition 
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for bistability, we start with the classical equations with detuning, A, included, 



— = e- {-f + iA)a + igAa{(3 + (3*), 

^ = -zUmP - lm(P ^9A\^. (19) 

These have the steady-state solutions, 



Pss PsS I'^SS] 1 



2 

2 



e 



which immediately leads to the cubic equation in J = la^sp, 



(20) 



/3 - — P + (72 + A^)I - = 0. 21) 

We find the condition for bistability by differentiating this expression with respect to /, 
which gives 

2 Sg'^A'^A , 2 .2^ 
- — / + (7^ + A^) = 0. (22) 

The condition for bistabihty is that this quadratic equation has two positive real roots. The 
roots are written as 

giving the inequality 

A±^V^2_3^2>o. (24) 

A necessary, but not sufficient condition is that A be positive, as / must be positive. This 
immediately contradicts the condition given by Mancini and Tombesi jsj, |A| > v^7, which 
allows for negative intensities. As J must also be real, we find the condition for bistability 
as A > 737. 

VI. STOCHASTIC RESULTS 
A. Initial conditions 

To numerically integrate Eq. we make use of the fact that the Ito and Stratonovich 
forms are identical so that we may use a standard three-step predictor-corrector method. 
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The convergence of the algorithm was checked by comparison with a four-step method, 
and also by halving the time-step. In all the quantities shown, the sampling errors are 
comparable to the thickness of the plotted lines. 

We will use the published experimental parameters of Ref. j^, and make comparisons 
with the theoretical predictions reported elsewhere. The oscillating mirror is considered 
as being perfectly reflecting, with a mass of m = 10~^ kg, a mechanical quality factor 
of Q = 4 X 10^ at 4.2 K, decreasing to 2.25 x 10^ at 70 K, and a resonance frequency 
uJm/2Tx = 26 kHz. The damping rate of the mirror is 7^ = -^uJm/Q = 0.0363 s^^ at 4.2 K. 
We consider a cavity length of L = 1 cm, with cavity finesse = 15 x 10^, which gives 
7 = 7rc/2jFL = 3.14 X 10^ s~^. We consider an optical wavelength of A = 1064 nm, which 
gives ujq = 1.77 x 10^^ s~^, and a coupling between the light and the mirror of g = ujq/L = 
1.77 X 10^^ m^^s^^. The optical pumping of the cavity is e = a/tPT^M); where V is the 
laser power in Watts. 

In stochastic integration of the equations which describe an intracavity optical system, 
the standard approach is to begin with the state inside the cavity as vacuum so that, 
with a continuous pump, the system enters the steady-state (or its limit cycle behaviour 
in the case of self-pulsing) after a few cavity lifetimes. In the present case the situation is 
somewhat different, as not only the electromagnetic field, but also the oscillating mirror, has 
to reach the steady-state. As the relaxation time of the mirror can be orders of magnitude 
larger than that of the intracavity field, and we need to average over a large number of 
trajectories to obtain reliable results, it is not practical to begin the integration with an 
arbitrary initial condition for the mirror. Naively beginning with /5(0) = jS'^iO) = 0, the 
ground state of the mirror, leads to extremely long lived transients, as this is far from the 
equilibrium state at finite temperature. To give some idea, even at 4.2 K, which is perhaps 
the lowest easily achievable temperature, the average number of mirror quanta becomes 
\P\ = ksT /hum = 3.36 X 10^. For purposes of comparison, we will therefore choose the 
initial mirror state in two different ways and integrate the equations without any pumping 
of the cavity. Firstly, as a (real) coherent state, which has uncertainties in position and 
momentum at the SQL, with the P-function 

P{(3) = 5{l3-^k^fjh^), (25) 
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and, secondly from the thermal distribution 

P(/5) = l-e-\P\Vn^ (26) 
vrn 

where n = ksT /hujm- Note that, at the beginning of each trajectory, (3 = {f3~^y and the 
phase is completely random for the thermal distribution. We stress here that the variance 
in the number of mirror phonons for a thermal state is V{n) =n'^ + n, which is very much 
larger than n. Using the Planck distribution, 

rl= (e^'^'"/'=^^-l)"\ (27) 

we find that to achieve n = 1, we would need T = 1.8 /iK, and even then the variance would 
be 2, or twice the mean value. 

It is important to note here that the number of phonons does not enter into the equations, 
but rather the quadratures and 1^. In a linearised approach using our equations, it is 
the uncertainties in these which are important. We can easily calculate these for a thermal 
state of the mirror with an unpumped cavity. A simple integration gives 

/ U nn\3/2 

ViX,) = ViY,) = 1 + 2^ ( ^ J , (28) 

equal to 2.2 x 10^° at T = 4.2 K for our system. This is in stark contrast to a coherent state 
of the mirror, sometimes used to facilitate the mathematics of a linearised analysis, and for 
which V{Xb) = V{Yb) = 1. 

We have calculated the stochastic results for the means and standard deviations of the 
mirror position, without any optical pumping and for initial thermal and coherent states 
of the mirror at temperatures of 4.2 K and 70 K. In Fig. ^ "we show the stochastic results 
for the position, x, of the mirror, with initial coherent states at these two temperatures. 
Considering only these mean values could give the erroneous impression that the mirror 
is in a non-stationary steady-state, which we can immediately see is not the case when 
we look at Fig. |21 which shows the standard deviations, cr{x), for the same parameters. 
Although we have shown the standard deviations here, the variances for the initial coherent 
state continue to increase linearly for more than twice the time shown, which was as far 
as we continued the integration. In contrast, for an initial thermal state, the mean value 
of the mirror displacement is, by definition, zero as can be seen from the equation for the 
P-function fj26|) . We note here that our stochastic results over more than 2 x 10^ trajectories 
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FIG. 1: Mean values for x at T = 4.2 K and T = 70 K for an initial coherent state of the mirror 
and no optical pumping. These results are the averages of 7.6 x 10'^ and 2.7 x 10^ trajectories, 
respectively. Note that, unless otherwise stated, the values plotted here and in subsequent graphs 
are dimensionless. 

still showed oscillations of the order of 10~^^ m, but that we are confident that this small, 
but non-zero, value is due to the difficulty of sampling the distribution with a finite number 
of trajectories. By comparison with the coherent state values, for an initial thermal state 
a{x) 1 X 10^^^ m at T = 4.2 K, and is almost constant, indicating that this is a good 
choice of initial condition. This value agrees well with the expression given in Ref. 1^ for 
the thermal noise, cr{x) = \/k^TJrmJ^, which gives a value of 1.47 x 10"^'^ m. Note that, 
over the time scales shown in Fig. ^ we do not see any decay in the oscillations towards the 
thermal state values, as this would be expected to happen on a time scale of l/7m, which is 



approximately 50 s for the parameters used here. In fact, although an initial co. 
of the mirror has been used in theoretical analyses (see, for example, Bose et al. 



lerent state 
) , it is not 



at all obvious how this particular state may be constructed experimentally. A thermal state 
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FIG. 2: Standard deviations for x at T = 4.2 K and T = 70 K for an initial coherent state of the 
mirror and no optical pumping. Note that these quantities were still increasing at twice the time 
shown here and would be expected to eventually attain the thermal values. 

arises naturally, and will be equal to a coherent state for T = K, but absolute zero cannot 
be reached experimentally. In optics, a coherent state can be described theoretically as a 
displacement of the vacuum by the displacement operator D{a) = exp(aa^ — a*a), which we 
can see has some relation to the optical pumping term of the Hamiltonian, Hpump- Therefore 
an ideal empty cavity with this pumping term will naturally develop an intracavity coherent 
state. We are not aware of any similar candidate for the mirror, even if it could begin in the 
T = vacuum state. Therefore we will use an initial thermal state in our investigations. 

When we examine the stochastic results for the intracavity field intensity, for an input 
power V = 5 mW, we find that the field exhibits a self-pulsing behaviour at approximately 
the resonance frequency of the mirror, as previously predicted However, the oscillations 
are of very small relative amplitude, at approximately 0.2% of the average mean intensity. 
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With increasing input power, the oscillations become larger until, for a power of lOOmW, 
for example, they are more than half the maximum intensity. At the lower power, the 
mean motion of the mirror is an oscillation between and 1.2 x lO^^^m, while at the higher 
power oscillates between —1 and 3 x lO^^^m. Even though these displacements are truly 
microscopic, they have a noticeable effect on the mean intensity, which should be easily 
detectable experimentally. Interestingly enough, these results are almost identical to what 
we find by numerical integration of the classical equations (fT^ . although these can tell 
us nothing about the quantum correlations which we wish to investigate. Among these 
quantum correlations are the variances of the intracavity field and the Fano factor, defined 
as F{Na) = V{Na)/Na- These results, which we averaged over 6.71 x 10^ trajectories, are 
shown in Fig. |21 For a coherent state, all three values are 1, which would be zero on the 
logarithmic vertical scale used here. As all three correlations are greater than or equal to 
one, we do not see any squeezing of the field in the time domain, but do see excess noise 
in all three quantities. As is common with Kerr media, there is more excess noise in the Ya 
quadrature than in either the Xa quadrature or in the intensity. 

B. Position measurements 

Much of the theoretical and experimental interest in this system has been in indirectly 
measuring small displacements or small forces which act on the oscillating mirror, by means 
of measurements on the optical field. Generally theoretical results are presented in terms of 
output spectra which allow for an inferred value of the mirror position at various frequencies. 
These spectra are simple to calculate in a linearised analysis which treats the system as an 
Ornstein-Uhlenbeck process Q], but we do not consider linearisation valid here, for the 
reasons we have stated. To calculate spectra from the results of stochastic integration is 
also possible in many cases, but here is made difficult by the stiffness of our equations, where 
the field and the mirror oscillate on vastly different time scales. The length of time needed to 
integrate a large enough number of trajectories over a sufficient time to give reliable results 
upon Fourier transformation is prohibitive. Hence we will show results which were obtained 
in the time domain. 

The first results we show for the mirror, in Fig. IH are for the uncertainties in the mirror 
position and momentum, defined as <y{x) = A^V{Xb) and a{p) = \B\^V(Yb). The interac- 
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FIG. 3: The variances of the intracavity electromagnetic field for a laser power of 5 mW. The solid 
line is V{Xa), the dash-dotted line is V{Ya), and the dotted line is the Fano factor. 

tion with the field has not noticeably changed these quantities from the values we found via 
stochastic integration in the unpumped thermal state case, but has given the momentum 
a mean value which oscillates between approximately ±1 x 10~^^ kgms~^, whereas it was 
stable at zero in the unpumped case. It is readily seen that the fluctuations in momentum 
are several orders of magnitude larger than the mean value. We also note here that the un- 
certainties shown are certainly underestimated, due to the difficulties of simulating the large 
(3 values with small probabilities in the initial conditions used. This result also shows that, 
when the likely thermal state of the mirror is taken into account, the back-action noise of the 
field on the mirror has very httle effect at the laser power considered here. At an increased 
power of 100 mW the standard deviations in both position and momentum both increase 
due to this back-action noise, reaching oscillatory values around means of ~ 1 x 10~^^, while 
the mean momentum oscillates between ±3 x 10~^^ kgms~^. 

We next investigate the degree of correlation between the position of the mirror, the 
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FIG. 4: The standard deviations in the mirror position and momentum for a laser power of 5 mW. 
The standard S.I. units for position and momentum are used. Note that a minimum uncertainty 
state of the mirror has a{x) = 5.6824 x 10^-'^'^ m and a{p) = 9.283 x 10^^*^ kgms~-^. 

intracavity optical intensity, and the Xa and Ya quadratures. This correlation function 
between two quantities w and z is defined as 

(wz) — {w){z) 



C{wz) 



(29) 



^V{w)V{z) ' 

where a perfect correlation gives a value of 1, a perfect anticorrelation gives a value of —1 and 
zero signifies no correlation at all. In Fig. El we see that, as predicted by Eq. [THl the strongest 
correlation is between Ya and the mirror position, these two being almost perfectly correlated. 
The correlation functions C{xXa) and C{xNa) oscillate at the frequency between a value 
of almost minus one and a value of approximately zero, ranging from being almost perfectly 
anticorrelated to almost perfectly uncorrelated. This behaviour is quite different from that 
shown in a linearised fluctuation analysis, where we found that all three of these functions 
showed anticorrelation around the mirror frequency and were zero at other frequencies. The 
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FIG. 5: The intracavity correlation functions, C{xXa) (solid line), C{xYa) (dash-dotted line), and 
C{xNa) (dashed line), for an input power of 5 mW. 

difference between our fully quantum nonlinear results and the usual linearised predictions is 
dramatically demonstrated when we consider ICfXaY;,)!^, used to analyse a possible quantum 
nondemolition measurement by Jacobs et al. [6j. Their prediction (see their Eq. 31) is for 
a maximum spectral value of 0.9757 for the parameters we use in this work, whereas our 
stochastic prediction in the time domain oscillates between and 6 x 10~^, showing almost 
no correlation at all. We think that it is unlikely that this difference can be explained by 
the difference between a spectral measurement and a time domain measurement, but that 
it is due to the inappropriateness of the linearisation procedure for this system. 

Another way to infer the mirror position is by linear estimation following measurements of 
the optical field. We follow a method proposed by Reid [3| in the context of a demonstration 
of the Einstein- Podolsky- Rosen paradox, and also outlined in Dechoum et al. js^. We 
assume that a measurement of the Ya quadrature allows for a linear estimate of X^, X^^* = 
cYa + d. This is consistent with the expansion of Ya given above, in Eq. ^| After optimising 
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for d, the RMS error in this estimate is given by 

= ((X, - cYa)^) - {X, - cYa)\ (30) 
which we may minimise as a function of c, finding 

v'(v;) 

We may then write 

V^-\X,) = V{X,) - i^^^^^, (32) 
a quantity which we may calculate via stochastic integration. The inferred uncertainty in a 



measurement of x will then be o"*"'^(x) = A\J V^^^ {Xb) , which will be equal to the standard 
quantum limit when the mirror is inferred to be in a minimum uncertainty state. We have 
also calculated the inferred position uncertainty using measurements of Xa and Na in the 
same manner. As shown in Fig. IHl and expected from the previous discussion, using the Ya 
quadrature gives better results than using either Xa or Na, which give oscillatory inferences. 
As it is, all of these give inferred uncertainties well above the SQL oi A = 5.68 x 10~^* m 
for the parameters used. We note here that the actual calculated value of cr(x), which 
oscillates between 1.025 x 10^^^ m and 1.03 x 10~^^ m, is greater than the value inferred 
through measurement of the Ya quadrature, which has a steady-state value of approximately 
5 X 10^^^ m. This is also the case with the inferred measurements of Ref. js^, and is due to 
the almost perfect correlation between the position of the mirror and the Ya quadrature. 



VII. NONCLASSICAL STATE PREPARATION 



The pendular cavity and variations on the theme have been proposed as useful devices 
for the preparation of nonclassical states of the cavity field and also of the mirror or mirrors. 
In this section we will give a brief review of some of these proposals, in chronological order, 
and explain why our results lead us to believe that they may not be as practical as the 
original authors suggest. 

Bose et al. Q suggest that this system may be used to prepare multi-component 
Schrodinger cat states of the field, near number states, and entangled states when two 
or more modes interact with the mirror, as well as Schrodinger cat-like states of the mirror 
prepared via quadrature measurements of the field. Beginning with a simple Hamiltonian 
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FIG. 6: The inferred uncertainties in the mirror position, A^jV^"^^ {X\^^ calculated as in Eg. 1321 
The solid line is the estimate using 1^, the dash-dotted line uses Xa and the dotted line uses N^- 

without any pumping or dissipation terms, they develop a time evolution operator which 
is seen to have a Kerr-like term, which leads them to believe that this system may exhibit 
some of the nonclassical features of a Kerr interaction. They then assume that both the 
cavity field and the mirror are initially in coherent states. While this simplifies the math- 
ematics, and is reasonably accurate for the cavity field, our calculations have shown that 
a coherent state is not a reasonable initial condition for the mirror. As the nonclassical 
features predicted depend on this initial condition, which at one stage of the paper is taken 
to be a vacuum state (for simplicity), and interaction of the mirror with the environment is 
considered to be equivalent to that of a cavity field with a zero temperature reservoir, our 
results lead us to believe that these predicted nonclassical states will be, at the very least, 
extremely difficult to observe. 

In a subsequent article the same authors propose that the pendular cavity may be 
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used to probe the decoherence of a macroscopic object, namely the mirror. In this work 
the authors accept that the mirror will begin in a thermal state and note that a mixture 
of Schrodinger cat states of the mirror can therefore be produced by the interaction with 
the field. It is assumed that a superposition of the Fock states |0) and |?7,) can be prepared 
inside the cavity, with the most simply prepared value of n = 1. Using the P-representation 
expansion of the density matrix in terms of coherent states, an expression is found for the 
time development of this density matrix. The field can then be measured via the interaction 
with a ground state two-level atom which is passed through the cavity. The probability of 
the atom being excited can be related to the decoherence rate of the spatially separated 
superposed coherent states of the mirror's motion. Examining this scheme with our pa- 
rameters, we find that, with n = 1, the spatial separation between two of the superposed 
coherent states has a maximum of Axmax = 1.4 x 10~^^ m, which is less than the thermal 
de Broglie wavelength, A^^ = 2.19 x 10^^^ m, so that it becomes difficult to claim we have 
a spatially separated superposition. When we take into account that the decomposition of 
the thermal state will include a huge number of coherent states, all with different phases, it 
seems that the practical realisation of this scheme may be more demanding technologically 
than the authors had supposed. 

Another proposal uses radiation pressure to entangle two macroscopic mirrors [l^, but 
relies on Langevin equations which are linearised around their classical steady-state solu- 
tions. As we have shown above, this process is of rather doubtful validity when the mirrors 
are coupled to thermal reservoirs. A further idea is to produce an Einstein-Podolsky-Rosen 
(EPR) state in the position and momentum of two spatially separated oscillating mirrors 
using the output of a nondegenerate optical parametric oscillator (OPO) [37J. We note here 
that, although the authors call the oscillator an optical parametric amplifier, the nonlinear 
crystal is inside a pumped cavity, so we will follow the usual terminology. (See their Fig. 2) 
Unlike most other treatments, this work uses an effective linear coupling between the light 
and the mirror. Like most of the others, the authors linearise Langevin equations around 
their steady state solutions. The OPO is treated via two-mode equations which do not 
describe the normal well-known threshold behaviour of such a device at all, and lead to the 
prediction of an entangled state of two combined quadratures, which is said to demonstrate 
an EPR correlation. The two output fields of the OPO are used to drive the two mirrors, 
which drives them into an EPR state of position and momentum. The mechanical damping 
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of the mirrors introduces a noise term Vb = 1 + 2nT, where Ht = coth{huJm/2kBT) is the 
mean thermal phonon number. This noise term seems to be the quadrature variance assum- 
ing a number state In-p), of the mirror. For the parameters we use, this expression gives a 
value ofVb = 6.7 x 10^, several orders of magnitude lower than our estimate in Eq. [2H1 which 
gives Vb = 2.2 x 10^°. Due to the issues we have raised here, we are lead to believe that the 
EPR state of the two mirrors may not be as easily demonstrated as indicated by the authors. 
Mancini et al. have also put forward a proposal to entangle two movable mirrors which 
form part of a four-mirror ring cavity. This work has the advantage of a description for the 
Brownian motion that is consistent with quantum mechanics at all temperatures, but the 
results presented are also obtained following a linearisation process, which we have shown 
to not be valid for reasonable temperatures. The final proposal we will consider here is by 
Marshall et al. , which treats the creation of superposition states of a macroscopic mirror 
{k. 10^^ atoms) via the interaction with a single photon. This work again uses a Hamiltonian 
approach without dissipation and assumes that the mirror can be prepared in its ground 
state, |0), which then allows for analytical solutions. The inclusion of decoherence and finite 

n 

temperature follow the approach of Ref. and we therefore have the same doubts about 
the physical viability of this proposal. 



VIII. COOLING BY FEEDBACK 



It may be thought that to reduce the thermal noise of the mirror it is sufficient to cool 
to a lower temperature and impressive reductions in thermal noise, of the order of 10^, have 
been achieved 10, 3^. However, it seems that the method of feedback cooling leads to a 
thermal equilibrium of the mirror at a lower temperature. For our system, using Eq. |^ 
a reduction by this factor leads to V{Xb) = yiXb) ~ 10^, which is still far above the 
coherent state level where linearisation of the equations may be expected to work. Ref. 
suggests that a mirror may be cooled to as low as 2 mK by dilution refrigeration, which 
would give us V{Xb) = V{Yb) ~ 2.27 x 10^, which, if we assume that feedback cooling from 
this temperature is as efficient as at room temperature, would allow the variances to be 
further reduced to something of the order of 10^. This is now actually smaller than the 
mean number of quanta, equal to 1.6 x 10^, although the variance in the number will 
now be of the order of 10^, so that linearisation will still not be reliable. 
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Standard cooling by feedback depends on the interaction of the electromagnetic field with 
the mirror phonons, which is proportional to a)aXb, which will not drive the mirror toward 
a coherent state. A coupling term proportional to o)h — may be expected to do this, 
while squeezing of the mirror position would conceivably be possible with a coupling of the 
type a)"^}) — a^6^ or h'^'^a — P'a). It is interesting to note that an effective coupling of the 
type required to drive the mirror toward a coherent state was used by Zhang et al. jsTj, 
although the physical parameters used to develop the coupling Hamiltonian are far from 
those we have used in our analysis. As an example in that work, uOm ^ 7, whereas we have 
used the values uom = 1.63 x 10^ s~^ and 7 = 3.14 x 10^ s"^. Zhang et al. also consider 
that the mirror damping is of the arnplitude form, which is not what we find using the 
Diosi master equation. Mancini et al. |4l| have actually presented a feedback scheme based 
on an effective Hamiltonian which couples a light quadrature with Xf,, later proposing that 
this method, which they call stochastic cooling, could be used to beat the SQL by achieving 
steady-state position squeezing of the mirror j^lf. The effective coupling used again depends 
on the linearisation of the equations of motion. We feel that whether feedback can be used 
to cool a macroscopic oscillator towards, or even beyond, the SQL remains an open question 
and subject to further research. 



IX. CONCLUSION 



In conclusion, we have presented a fully quantum analysis of the pumped pendular cavity 
which does not depend on linearisation of the equations of motion. We have shown that 
the thermal noise of the mirror is overwhelming for typical temperatures and experimental 
parameters. This means that the linearisation procedure commonly used is of doubtful 
validity. Unless ways can be found to change the quantum state of the mirror in respect 
to the thermal excitations, we expect that the nonclassical states predicted in a number of 
theoretical analyses will not be able to be demonstrated experimentally. Feedback cooling 
techniques, while able to lower the effective temperature to an impressive degree, serve 
merely to produce another thermal state of the mirror at a lower temperature. The noise is 
still overwhelming if we wish to reach the SQL, or even beat it, which would be necessary 
for the detection of gravity waves. This is also the case if we wish to observe some of the 
quantum superpositions and entanglement which have been theoretically predicted. 
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